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^ ABSTRACT 

Aims. We study the evolution of the luminosity function (LF) of type-1 and type-2 Active Galactic Nuclei (AGN) in the mid-infrared, and then 
derive the contribution of the AGN to the Cosmic InfraRed Background (CIRB) and the expected source counts to be observed by Spitzer at 
24yum. 

Methods. We used a sample of type-1 and type-2 AGN selected at 15 /jm (ISO) and 12yum (IRAS), and classified on the basis of their optical 
^ spectra. Local spectral templates of type-1 and type-2 AGN have been used to derive the intrinsic 15 fim luminosities. We adopted an evolving 

smooth two-power law shape of the LF, whose parameters have been derived using an un-binned maximum likelihood method. 
CZJ Results. We find that the LF of type-1 AGN is compatible with a pure luminosity evolution (L(z) = L(0)(l + zf'-) model where ~2.9. A 
small flattening of the faint (L15 < LJj) slope of the LF with increasing redshift is favoured by the data. A similar evolutionary scenario is 
^ ^ found for the type-2 population with a rate ranging from ~1.8 to 2.6, depending significantly on the adopted mid-infrared spectral energy 
• ■ distribution. Also for type-2 AGN a flattening of the LF with increasing redshift is suggested by the data, possibly caused by the loss of a 
^\ ' fraction of type-2 AGN hidden within the optically classified starburst and normal galaxies. The type-1 AGN contribution to the CIRB at 15 yum 
^ , is (4.2-12. l)x 10"" Wm^^sr"', while the type-2 AGN contribution is (5.5-1 1.0)x 10"" Wm"^ sr"'. We expect that Spteer will observe, down 
■ - - I to a flux limit of 524^m=0.01 mjy, a density of ~1200 deg"- type-1 and ~1000 deg"- type-2 optically classified AGN. 

Conclusions. AGN evolve in the mid-infrared with a rate similar to the ones found in the optical and X-rays bands. The derived total contri- 
bution of the AGN to the CIRB (4-10%) and Spitzer counts should be considered as lower limits, because of a possible loss of type-2 sources 
caused by the optical classification. 
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1. Introduction 

Active Galactic Nucleus (AGN hereafter) are expected to have 
played an important role in the formation and evolution of the 
galaxies in the Universe. An example is the observation of the 
correlation between the mass of the central black holes (Mbh) 
and the mass of the bulges (Magorrian et al. 1998) or the veloc- 
ity dispersion of gas and stars in the bulges of galaxies (Mbh-o' 
relation; Ferrarese & Merritt 2000; Tremaine et al. 2002). Thus, 
the measure of the history of the density of AGN as a function 

* Data presented here are based on observations by the Infrared 
Space Observatory (ISO). ISO is an ESA project with instruments 
funded by ESA member states (especially the PI countries: France, 
Germany, the Netherlands and the United Kingdom) and with the par- 
ticipation of IS AS and NASA. 



of the luminosity (the luminosity function, LF hereafter) can 
provide fundamental clues to explain the present day universe 
(e.g. Balland et al. 2003; Menci et al. 2004; Granato et al. 2004; 
Di Matteo et al. 2005). 

The LF will not only provide informations on the demo- 
graphics of AGN, but will also place constraints on the phys- 
ical model of AGN, the origin and accretion history into su- 
permassive black hole and the formation of structures in the 
early universe. Moreover, the measure of the AGN LF over the 
whole spectral range will provide information on the relative 
importance of accretion power in the overall energy budget of 
the universe. 

AGN have been historically classified into two groups 
(type-1 and type-2), depending on the presence or absence in 
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their optical spectra of broad emission lines (FWHM>1200 
km s '). The unified model for AGN assumes that the two types 
of sources are intrinsically the same (Antonucci et al. 1993), 
and that the observed differences in the optical spectra are ex- 
plained by an orientation effect. The presence of obscuring ma- 
terial (torus or warped disc) around the central energy source 
and its orientation respect to the observer could prevent a di- 
rect view of the central region around the nuclei responsible 
for the broad line emission observed in the optical band. AGN 
are then classified in two groups: the un-obscured (type-1) and 
obscured (type-2) sources depending whether the line-of-sight 
intersects or not the obscuring material. However, nowdays it 
seems clear that the orientation-dependent model is probably 
just a Oth-order approximation to the true nature of AGN (e.g. 
Barger et al. 2005; La Franca et al. 2005). 

The major efforts to understand the evolution of AGN have 
been historically concentrated at optical wavelengths and in the 
last years have produced the largest compilations of AGN from 
the Two Degree Field (2dF; Boyle et al. 2000; Croom et al. 
2004) and the Sloan Digital Sky Survey (SDSS; York et al. 
2000; Hao et al. 2005). However the optical bands have proved 
to be very inefficient in the selection of obscured sources and 
the assembly of unbiased samples of type-2 objects is diffi- 
cult. In the last years many new windows have been opened at 
various wavelength bands for the observation of high redshift 
galaxies and AGN: sub-mm (SCUBA), Infrared (IRAS, ISO and 
Spitzer) and the X-ray (ROSAT, XMM and Chandra). Thanks 
to their smaller dependence on dust obscuration if compared 
to the optical surveys, the Infrared and the hard X-ray wave- 
lengths have proved to be much more efficient in the detection 
of type-2 sources. 

X-ray observations have been so far consistent with 
population synthesis models based on the Oth-order unified 
AGN scheme (e.g. Comastri et al. 1995) or its modifications 
(PompiHo et al. 2000; GilH et al. 2001; Ueda et al. 2003; La 
Franca et al. 2005). The hard spectrum of the X-ray back- 
ground is explained by a mixture of absorbed and unabsorbed 
AGN, evolving with cosmic time. According to these models, 
most AGN spectra should be heavily obscured as the light pro- 
duced by accretion is absorbed by gas and dust. 

For this reason an important mid-infrared (4-40 jjm) ther- 
mal emission is expected from reprocessed radiation by gas 
and dust grains directly heated by the central black hole 
(e.g. Granato et al. 1997; Oliva et al. 1999; Nenkova et al. 
2002). AGN could be in this case important contributors to 
the Cosmic Infrared Background (CIRB). Indeed, IRAS has 
observed strong mid-IR emission in all local (z < 0.1) AGN 
(Miley et al. 1985; de Grijp et al. 1985; Neugebauer et al. 1986; 
Sanders et al. 1989). Complete and largely unbiased samples of 
AGN have been produced at 12 jjm (Rush, Malkan & Spinoglio 
1993, RMS hereafter) and 25 (Shupe et al. 1998) based on 
IRAS observations. These samples provide an important under- 
standing of the infrared emission from local galaxies and a firm 
base to compare the properties of galaxies in the local universe 
with the high redshift populations uncovered by ISO (z < 1) 
and most recently by Spitzer (z < 2). 

Statistical studies of AGN in the mid-IR have been based on 
large but local samples from the IRAS observations (e.g. RMS, 



Shupe et al. 1998) and deep but small ISO fields in the Hubble 
Deep Field North (HDFN; Aussel et al. 1999a, 1999b), South 
(HDFS; Oliver et al. 2002) and the Canada-France Redshift 
Survey (CFRS; Flores et al. 1999) that provided only a handful 
of objects. 

Therefore, due to the small number of high redshift objects 
available, the shape and evolution of the mid-IR LF of these ob- 
jects at high redshift was largely unknown. Only recently, large 
and highly complete \5 spectroscopic samples, based on 
the European Large Area ISO Survey (ELAIS hereafter; Oliver 
et al. 2000), have been released (Rowan-Robinson et al. 2004; 
La Franca et al. 2004; FLF04 hereafter). Based on the analysis 
of the ISO observations in the ELAIS-Sl field, Matute et al. 
(2002, M02 hereafter) derived the first estimate of the type-1 
LF in the mid-IR (15yum) finding an evolution similar to the 
ones observed in the optical (Croom et al. 2004) and X-rays 
(e.g. Hasinger et al. 2005; La Franca et al. 2005). The evo- 
lution of the LF for normal and starburst galaxies detected at 
15jum in the southern ELAIS fields (SI & S2) has been derived 
and discussed by Pozzi et al. (2004). 

In this paper we investigate the separate evolution of type- 
1 and type-2 AGN in the mid-IR and their contribution to the 
CIRB using: a) the local IR population uncovered by IRAS, 
b) deep observations from ISO fields and, c) the most recent 
spectroscopic identifications of \5pm sources in the southern 
fields of the ELAIS Survey (Pozzi et al. 2003; FLF04)i. The 
AGN sample used in our work is described in section 2. Section 
3 introduces the general method used to derive the parameters 
of the luminosity function and presents the results. In section 4 
we discuss our results and compare them with the LFs derived 
in the mid-IR by previous analysis of mid-IR observations. The 
contribution to the CIRB and the predicted counts in the mid-IR 
Spitzer band at 24 pva are computed and discussed in sections 
5 and 6. Section 7 summarises the main results. 

A cosmology with Hq - 75 km s"' Mpc"' and Q.x - 0.7 , 

= 0.3 was adopted for the work presented in this paper. 

2. The sample and source classification 

The mid-IR selected AGN sample used in this work has been 
extracted from: (/) the 15 yum ELAIS fields SI and S2 (Lari et 
al. 2001; Pozzi et al. 2003), (//) the 15 pm deep ISO surveys in 
the HDFN (Aussel et al. 1999a, 1999b), HDFS (Oliver et al. 
2002) and CFRS (Flores et al. 1999), (Hi) the local IRAS 12 pm 
sample of RMS. 

The 15jum catalogue in the ELAIS SI field has been re- 
leased by Lari et al. (2001). It covers an area of ~4 deg^ cen- 
tred at Q-(2000) = 00''34'"44.4', 5(2000) = -43''28'12" and 
includes 462 mid-IR sources down to a flux limit of 0.5 mJy. 
Mid-IR source counts based on this catalogue have been pre- 
sented and discussed by Gruppioni et al. (2002). 

The analysis presented here is restricted to the more reli- 
able subsample of 406 sources as described by FLF04. About 
80% of these sources have been optically identified on CCD 
exposures down to R~23, while spectral classification has been 

' Data and related papers about the ELAIS southern survey are 
available at: http://www.fis.uniroma3.it/~ELAIS-S_ 
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Fig. 1. Distribution of ELAIS-Sl type-1 and -2 AGN in the z. - /^-magnitude plane. 
Best fit linear solution (±lo") in the log(z)-R space are shown as dotted lines for type-1 
AGN and as dashed lines for type-2 AGN. The -magnitude and redshift distributions are 
shown in the left and bottom panels respectively. 

obtained for 90% of the optically identified sample. As dis- 
cussed by FLF04, due to a different mid-IR flux limit coverage 
of the ELAIS-Sl field, the total area has been divided into two 
regions: the central and deepest region of SI (Sl-5) reaching 
mid-IR fluxes (5 15 hereafter) of 0.5 mJy, and the outer region 
(SI -rest) with a 0.9 mJy flux limit. The Sl-5 area is spectro- 
scopically complete at the 97% level down to /?=21.6, while 
Sl-rest completeness reaches the 98% level down to /?=20.5. 
In total 116 sources (29% of the total mid-IR sample) do not 
have a spectroscopic identification due to incompleteness of 
the follow-up or to the lack of optical counterpart brighter than 
/?=23. A detailed description of the optical identification, spec- 
troscopic classification, size and completeness function of the 
different areas used in the ELAIS-Sl sample are presented and 
discussed by FLF04. 

The S2 field is a smaller and deeper area centred at 
Qr(2000) = 05''02'"24.5*, 5(2000) = -BO^Se'OO" covering ~ 
0. 12 deg^. This region includes 43 sources with SfN > 5 down to 
S 15=0.4 mJy (Pozzi et al. 2003; Rowan-Robinson et al. 2004). 
Photometry in the whole field is available in the U, B, I and K' 
bands down to 21.0, 24.5, 22.0 and 18.75 respectively (Vega 
magnitudes). 39 infrared sources have a counterpart in the / 
band. 22 of them have a spectroscopic identification, while 8 
sources are classified as stars since they are associated with 
bright (/<14) point-like sources in the optical catalogue. To 
avoid large incompleteness due to the optical follow-up, only 
sources with optical counterparts brighter than 7=20.6 were se- 
lected. Down to this optical limit, the sample considered in S2 
is then complete at the 95% spectroscopic level for sources with 
S 15 >0.7 mJy. 

Classification for AGN dominated sources in the ELAIS 
fields was based on their optical spectral signatures. Sources 
showing broad emission line profiles (rest frame FWIIM> 1200 
kms^^) were classified as type-1. Type-2 sources were se- 
lected following classic diagnostic diagrams (e.g. Tresse et 
al. 1996; Osterbrock 1989; Veilleux & Osterbrock 1987) 
that included one or more of the following Une ratios: 
[Nn]/Ha, SII/(Ha-i-[Nn]), OI/Ha, [Om]/I^ and [OII]/I^ 



when available, depending on the redshift of the source {e.g. 
log([OIII]/I^)>0.5 and log([NII]/Ha)>-0.4). 

In total, the southern ELAIS (S1+S2) AGN sample in- 
cludes 27 (25-H2) type-1 and 25 (23+2) type-2 AGN represent- 
ing ~16% (52/320) of the identified infrared sources and -24% 
(52/221) of the extragalactic population. Their distribution in 
the redshift-/?-magnitude space is shown in Fig. 1 . The plotted 
lines represent a linear fit (+1(t) to the data {dotted for type-1 
AGN and dashed for type-2). Type-1 sources are detected up 
to z ~3 and show a quasi redshift-independent and narrow op- 
tical magnitude distribution {R = 18.76 + 1.01). On the other 
hand, type-2 sources show a larger spread of optical magni- 
tudes (/?= 18.45 ± 1.62) and a rather steep dependence on the 
redshift up to z~l. 

Constraints on the fainter population are provided by deep 
15 jum observations in the HDFN (Aussel et al. 1999a, 1999b), 
HDFS (Ohver et al. 2002) and in the 1452-1-52 field of the 
CFRS (Flores et al. 1999). These fields have a flux limit about 
an order of magnitude deeper (5 15 - 0.1 mJy) than the ELAIS - 
S fields but with sky coverages 100 to 500 times smaller. The 
multiwavelength foUow-up in these fields has identified a small 
fraction of sources as AGN (both type-1 and type-2). 

In the HDFN, 41 15jum sources have been detected over 
an area of 21.5 arcmin^. The sample is complete down to 0.1 
mJy (Aussel et al. 1999a). 4 of these sources, with fluxes be- 
tween 0.1 and 0.45 mJy, are found to be AGN by Alexander et 
al. (2002), equally divided between the two classes (Alexander 
et al., private communication). In the HDFS we use the AGN 
identified by Franceschini et al. (2003) from a sample of 59 
ISOCAM sources with 5 15 > 0. 1 mJy selected over an area of 
19.6 arcmin^. They have detected 2 secure AGN, one type-1 
(Si5 = 0.288 mJy) and one type-2 (S15 = 0.518 mJy). The ob- 
servations in the CFRS 1452-1-52 field cover an area of ~100 
arcmin^. 41 sources with 5 15 > 0.35 mJy and with S/N>4 have 
been selected (Flores et al. 1999). The spectroscopic follow- 
up has identified one type-1 {S 15 = 0.459 mJy) and one type-2 
{S 15 = 1.653 mJy) source. In the above mentioned fields, a frac- 
tion of the type-2 sources were not classified according to their 
optical spectra but on the basis of the shape of their multiwave- 
length spectral energy distribution (SED). For this reason only 
the selected type-1 population has been used in our analysis. 

Finally, we have combined the ISOCAM data with the local 
IRAS sample of AGN from RMS. The RMS Catalogue is a high 
galactic latitude {\b\ > 25°), colour selected (>S'60/im/>S'i2/im) and 
flux limited {S n^im >0.22 Jy) sample extracted from the IRAS 
12 nm Faint Source Catalog, Version! (Moshiretal. 1991). To 
avoid completeness uncertainties in the mid-IR sample, only 
sources with S i2^m >300 mJy were selected (see RMS). In to- 
tal 41 type-1 and 50 type-2 sources are found within this flux 
range. The photometric and spectroscopic follow-up of the op- 
tical counterparts is considered to be ~ 100% complete for these 
sources and has been used in our analysis as representative of 
the mid-IR selected population of AGN in the local universe. 
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Table 1. Mean redshift and Luminosity (±lo-) of the mid-IR 
samples 



Sample 



< z> 



IRAS 

<logL,5 



<z> 



ISO 



Type-1 0.03 ±0.02 
Type-2 0.02+0.01 



43.87±0.87 1.23±0.71 45.22±0.78 
43.66±0.68 0.28+0.20 43.90+0.69 



(a) 



12 fim luminosities from RMS have been converted to 15 fim using the adopted SED. 
In particular the Circinus SED was adopted for type-2 AGN (see §3.1). 
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Fig. 2. Adopted mid-IR spectral energy distributions (SEDs) for type-2 sources {top 
panel) and the corresponding /^-corrections as a function of redshift computed for the 
IS0CAM-LW3 filter (bottom panel). 



3. The Evolution of AGN 

3.1. The Method 

In order to compare local sources with a higher redshift se- 
lected population and to probe changes with cosmic time it is 
necessary to compare their intrinsic physical properties. Our 
sample has been selected in the mid-IR and identified in the 
optical band, therefore rest-frame luminosities, in the form of 
vLy, at 15 fim (L15) and in the /?-band (L^) were computed us- 
ing well known SEDs for type-1 and 2 sources. 

The compilation of Elvis et al. (1994) of 47 QSO, in the 
range of 1-20 /.<m, was used as representative of the mid-IR 
emission for the type-1 population. This composite SED agrees 
very well with spectroscopic observations of type-1 sources 
performed by ISO in the mid-IR (Clavel et al. 2000; Spoon et 
al. 2002), showing a strong power-like continuum, very weak 
or no PAH emission bands and no evidence of the 10 //m sil- 
icate absorption feature. Optical luminosities for type-1 AGN 
were computed using the /?-band A:-correction by Natali et al. 
(1998), which takes into account the large effect of broad lines 
entering and leaving the passband. 



Unlike type-1 sources, the mid-IR SED of type-2 objects 
varies greatly. They range from starburst-like galaxies SEDs, as 
in the case of Circinus (Sturm et al. 2000), showing a weak con- 
tinuum, prominent emission from PAH molecules and a deep 
silicate absorption feature at ~10yum, to more power-like SEDs 
dominated by hot dust directly heated by the active nucleus, as 
in NGC 1068 (Sturm et al. 2000). As a consequence, the mid- 
IR SED of NGC 1068 and Circinus galaxies (Fig. 2) were used 
as representative of two extreme cases of obscured AGN in the 
mid-IR. In order to derive the optical ^-correction for type-2 
sources, an average spectrum was produced in the optical band 
(3500-7000A) using our best type-2 spectra (Fig. 5 in FLF04). 

In Figure 3 we show the distributions in the L15 - z plane of 
the total sample used in our analysis, while Table 1 summarises 
the mean redshifts and mid-IR luminosities, and corresponding 
Icr dispersion, measured for the local (IRAS) and high redshift 
(ISO) sources. 
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Fig, 3. Distribution in the L15 - z space of the total sample of mid-IR type-1 and 2 
AGN. The mid-IR luminosity for type-2 sources has been computed applying the Circinus 
^-correction (see text, §3.1). The dotted, dashed and dot-dashed lines represent the flux 
limits at 0.1, 1.0 and 300 mjy respectively. 



The adopted shape of the luminosity function (LF) is a 
smooth two-power law of the form: 



^flogLis (Z.i5/Z,;5)«fe) + {U,IL\^f 



(1) 



where we have introduced a dependence on z in order to 
allow a change in the slope (or) of the faint end of the luminosity 
function. We have assumed a luminosity evolution of the form: 



il5(0)(l 



•z)*^ for z<Zcut 
for z > Zcut 



(2) 



where Zcm is the redshift cut beyond which the luminosity func- 
tion is assumed to remain constant. 

The parameters for the luminosity function and the evolu- 
tion have been derived using an un-binned, maximum likeli- 
hood method (Marshall et al. 1983), as described by Matute 
et al. (2002). The spectro-photometric follow-up of the ELAIS 
southern fields was not deep enough to provide an identification 
and a classification for all the sources detected in the mid-IR 
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sample. A factor depending on redshift and mid-IR luminos- 
ity was therefore introduced in order to take into account the 
fraction of objects not identified, due to the lack of an optical 
counterpart or to incompleteness of the spectroscopic follow- 
up. We will refer to this factor as the optical - completeness 
factor, 0(z,L). It represents the probability that a source with 
a given redshift and mid-IR luminosity L15 has an apparent R- 
band magnitude within the spectroscopic limits of the samples, 
and was derived taking into account both the average Lis/Lr 
ratio of the sources and its observed natural spread. For any 
given mid-IR luminosity L15, this probability has been com- 
puted assuming a gaussian distribution of log (Ln/Lg) centred 
at its mean value and with a sigma equal to the observed natu- 
ral spread of the Lis/Z-s relation (see next sections for discus- 
sions on the adopted mean values and spreads of log (Lij/Ls)) 
^. Then, the function to be minimised can be written as: 

S = -2 g ln[(l)(z,-,L,)] + 2jJ ^(z,L)Q.(z,L)&(z,L)^dzdL (3) 

where Q(z,L) is the effective area covered by each subsample 

at 15 pm, L-L\s and the index ; runs over all the objects of 
a given class. For the ELAIS-Sl sources the area coverage is 
presented as a function of flux by FLF04. Area coverages for 
the deeper ISO fields (HDF-N, -S and CFRS) are described in 
the previous section. The area sampled by RMS is considered 
equal to for 5 12 < 300 mJy and equal to 22191.80 deg^ for 
5 12 > 300 mJy (see RMS for details). The representative SEDs 
were used to convert the RMS 12 pm fluxes into 15 pm fluxes. 

Confidence regions for each parameter have been obtained 
by computing AS {- A;^f^) at a number of values around the best 
parameter (Smin), while allowing the other parameters to float 
(see Lampton, Margon & Bowyer 1976). A 68% confidence 
region for any single parameter corresponds to A5 = 1 . The 
goodness of the fit has been verified using the bidimensional 
Kolmogorov-Smirnov test (2D-KS hereafter) in the L15-Z space 
(see Peacock 1983; Fasano & Franceschini 1987). The normal- 
isation factor O* is determined in such a way to reproduce the 
observed total number of sources (ISO + IRAS). 

3.2. The Evolution of type-1 AGN 

For type-1 sources, the optical completeness factor, 0(z, L15), 
has been estimated using the mean value log(Li5/LR)=0.23 
with a Icr dispersion of 0.25 (from the sample of Elvis et 
al. 1994). These values are in rough agreement with the de- 
rived value for the IRAS local sub-sample of type-1 AGN from 
Spinoglio et al. (1995; 0.05 with a Icr spread of 0.13). This ra- 
tio and its spread shows no significant dependence on L15 over 
more than 4 decades of mid-IR luminosity and we assumed 
them to be constant with redshift. 

In order to find the best fit solutions to the LF we have 
considered: 

- 2 different behaviours for the faint slope of the LF: a) a 
fixed value (a = a(0)) and, b) a dependence of or on z of the 

^ A similar approach was adopted in the previous estimates of the 
MIR type-1 AGN LF by Matute et al. (2002), and of the X-ray AGN 
LF by La Franca et al. (2002, 2005) 



form: a(z) = ai exp(-a2 * z) -1- £1^3. A justification for this 
form is given in the following paragraphs. 

- A fixed value for Zcut equal to 2.0 in our parameterisation of 
the LF. Our relatively small number of sources, especially 
at high z, does not allow to sample correctly the existence of 
this redshift cutoff. We have therefore assumed a behaviour 
in the mid-IR similar to what is already known for the evo- 
lution of AGN in the optical (Boyle et al. 2000, Croom et 
al. 2004) and X-ray (Miyaji et al. 2000, 2001; La Franca et 
al. 2005; Hasinger et al. 2005), suggesting values for Zcm in 
the range 1.7-2.5. 

- A universe volume between z=0 and z=4, while the 
space density was integrated in the luminosity interval 
logLi5=[42,47]. 

In total, 72 (27 ELAIS -1- 41 RMS -1- 4 750-Deep) type-1 

AGN were used to derive the LF. The best solutions found for 
each parameter (+lcr confidence levels) and the 2D-KS prob- 
abilities are presented in Table 2 in the rows labeled "A" and 
"B" for the fixed and variable faint slope models respectively. 

In Figure 4 we show the observed space density distribution 
and the best fit models in two redshift intervals: z = [0,0.2], 
where the IRAS sources dominate, and z = [0.2,2.2], mainly 
populated by ISO sources. The large interval at high redshifts 
(z = [0.2, 2.2]) was chosen only for representation purposes, 
to assure a significant number of observed sources in each lu- 
minosity bin. In order to correct for evolution within the red- 
shift intervals, the observed space density distribution is plot- 
ted following La Franca et al. (1997). The expected number 
of sources given by the model in each bin (j\)'"'"^<^'^) is com- 
puted and compared to the observed number of AGN in the 
bin {N'^"'"). The ratio Ndataij^modei f^j. g^^j^ i^jj^ ^^^^ j^^jjj. 

plied for the value of the LF at the corresponding central lu- 
minosity and redshift value of the bin. The plotted error bars 
correspond to the Icr Poisson distribution and were estimated 
following Gehrels (1986). Space density upper limits are given 
where sources are expected by the model but not observed. The 
Icr dispersion for the LF was computed according to the Icr 
uncertainties of its parameters, and represented as a light-grey 
shaded area in Fig. 4. 

Beside the natural statistical uncertainties due to the limited 
number of sources used in our analysis, our fitting technique is 
affected by possible errors due the assumed average value of 
the LisILr ratio used to correct for the spectroscopic incom- 
pleteness of the samples. In order to estimate these uncertain- 
ties we have recomputed the LF best fit solutions assuming two 
extreme cases for the log(Li5/LK) ratio. In particular, we have 
used as average value of log(Li5/LR) our best estimate (0.23) 
plus or minus the observed Icr dispersion (0.25). Even under 
these extreme assumptions, the parameters of the LF result to 
be within the errors quoted in Table 2. The dark-grey shaded 
areas in both panels of Fig. 4 show the corresponding Icr range 
of uncertainty introduced in the estimate of the space density 
of type-1 AGN. As it is clearly seen in the figure, the statistical 
errors due to the limited number of sources used in our analysis 
dominate over the uncertainties due to the assumptions on the 
average value of the log(Li5/LR) ratio. For example, the uncer- 
tainty introduced in the estimated density of type- 1 AGN by the 
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Table 2. Parameter values of the fit of the Luminosity Function 
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0.33 
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6.00 ** 


76—" 

"•'"-0.14 

0.00 ** 
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" Ljg corre.sponds to vLy at 1 5 /jm and is given in erg s ^ . 
The normalisation, log , is given in Mpc"^ mag"' . 

Contribution vly to the CIRB at 15/im in units of 10"" Wm"-sr"' . In parenthesis we report the % of the CIRB produced by the fit given a CIRB value of 2.7 nW m"^ sr"' . 
** The values of these parameters reached the physical limit imposed in our minimisation process. These constraints were imposed to the models in order to not produce unphysical 
solutions at high redshift due to the low statistics. Therefore, no errors are reported for these parameters. 
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Fig. 4. Type- 1 AGN local (z=O.I, continuous line) and high redshift fe=1.2, dash- 
dotted line) best fits to the 1 5 //m Luminosity Function assuming a redshift dependent faint 
slope (Fit "B", top panel) and a PLE model (Fit "A", bottom panel). Open circles represent 
the observed space density of the local population (dominated mainly by RMS sources) 
while filled circles represent the observed space density of the high redshift population 
(dominated by ELAIS sources). Poissonian errors at lo" confidence level are shown. The 
dark-grey shaded areas (quite narrow and thus barely visible) represent the uncertainty 
introduced by the assumptions used to estimate the optical incompleteness of the samples, 
while the Mght-grey shaded areas correspond to the la- errors of the best fit parameters of 
the LP (see §3.2), 



assumed Lis/Lr relation is ~10% at Lis ~ L'^s' while the statis- 
tical uncertainty due to the limited number of sources is larger 
than 20%. For these reasons, in order to compute the errors 
over the parameters of the LF and any derived quantity (e.g. 
integral counts and redshift distributions), we neglected the un- 
certainties related to the assumptions on the average value of 
the log(Li5/Ls) ratio. 

We observe that although the overall fit by the PLE model 
(model A in Table 2, Fig. 4 bottom panel) is in reasonable 
agreement with the observed data, it overpredicts by a small 



amount (~lo"), but systematically over 5 luminosity bins, the 

number of the low luminosity (Li5(z) < L*^) and high redshift 
(z - 1.2) type-1 AGN. This is the main reason why a depen- 
dence of the faint slope on the redshift was included. The intro- 
duction of this dependence translates into a luminosity depen- 
dent luminosity evolution (LDLE) for the faint part of the LF 
(model B in Table 2). In this case, a better agreement between 
the model and the data is obtained (Fig. 4, top panel) as shown 
by the increment of the 2D-KS probability (0.35 vs. 0. 11). 

The measured rate of luminosity evolution (kL) results to be 
rather independent from the adopted relation between a and z. 
The best fit values, ki ~ 2.7-2.9, are similar, within the errors, 
to those already found for these objects in the optical and X- 
ray wavelengths, where the evolution rate is between 2.5 and 
3.5 (La Franca et al. 1997; Boyle et al. 2000; Miyaji et al. 2000 
& 2001; Croom et al. 2004; La Franca et al. 2005; Hasinger et 
al. 2005). Therefore, according to this analysis, no difference 
is found in the evolution of mid-IR selected type-1 AGN in 
comparison to the ones selected in the optical and X-rays. 

The 15 /im integral counts and the observed redshift dis- 
tribution for type-1 AGN can now be compared with the pre- 
dictions derived from the best fit models. The results of this 
comparison are presented in the upper right and left panels of 
Figure 5. The bottom panels show the corresponding relative 
errors as derived from the Icr uncertainties on the best fit pa- 
rameters of the LF. 

The expected integral counts (Fig. 5, left) given by both 
models {dashed line for model "A" and continuous Une for 
model "B") provide a good representation of the observed 
data (shaded areas) down to ~2 mJy. At fluxes below ~2 mJy 
the variable slope model agrees better with the ELAIS data, 
while the PLE model is more representative of the fainter pop- 
ulation in the HDF and CFRS fields. At the faintest fluxes 
(515 = 0.1 mJy) the variable slope model underpredicts the 
observed counts by a factor ~2. The uncertainties on the pre- 
dicted counts at this flux (~20%) are not large enough to justify 
the observed discrepancy. 
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Fig. 5. Left ) Mid-IR integral counts for type-1 AGN. Shaded areas show the RMS and ELAIS observed counts, while symbols (triangle, star and diamonds) represent the /50-Deep 
surveys. The thick solid line gives the predicted counts from the variable slope fit model ("B"). Expectations from the fixed faint slope model ("A") are plotted as a dashed line. All 
errors in the observed distributions are quoted at the Icr level and given by Poisson statistics. The bottom panel show the lo" relative error on the predicted counts due to the uncertainties 
on the best fit parameters of the LF (see text). Line types have the same meaning as in the above panel. Right) Redshift distribution for type-1 AGN. The liistogram represents the total 
number of observed sotirces in each redsliift interval. Shaded regions outline the ISOCAM contribution, while the open histogram accotmts for the contribution from IRAS sources. 
The lines are the predictions from the best fit model ("B"): thick dash-dotted for IRAS, thick dashed for ISO population, while the continuous thick line gives their sum. The expected 
distributions for the "A" model are plotted as above but using thin lines. The bottom panel shows the relative errors of the predictions for the entire mid-IR population (IRAS+ISO) only. 



We also find an overall good agreement between the red- 
shift distribution of the observed data and the model predic- 
tions (Fig. 5, right). The largest discrepancy between the mod- 
els occurs at the lowest redshift bin (z - [0, 0.5]) where the 
fixed faint slope model (thin lines) overpredicts the number of 
sources observed in the ISOCAM fields by a factor ~3 (4 ob- 
served vs. ~12 expected), while the variable faint slope model 
provides a better representation of the observed distribution. 

In summary, although the variable faint slope model pro- 
vides a better overall fit to the LF and reproduces better 
the z-distribution, it underestimates the faint mid-IR counts 
(5 15 < 0.5 mJy) by a factor of ~4. As both models ("A" fixed 
& "B" variable faint slope) are statistically acceptable, we con- 
clude that given the presently available statistics the two mod- 
els are equivalent within the errors. 

M02 measured for the first time the evolution of type- 
1 AGN based on a preUminary mid-IR catalogue from the 
ELAIS-Sl field. The values found here for the fixed faint slope 
LF (model "A") are in agreement with what already found by 
M02 within ~lcr errors. The only exception is for the bright 
slope (J3), where the here derived value is flatter (2.29^^ jg vs. 
2.89^g 2g). This difference can be understood if we take into 
account the large uncertainties in the brighter luminosity inter- 
vals (in both M02 and our sample) due to the small number of 
sources. While the 2D-KS test gives a 11% probability that the 
our sample of AGN is well represented by a non-evolving faint 
slope (model "A"), this probability decreases to 5% if we use 
the parameters values reported by M02 (Table 1 in M02). 

3.3. The Evolution of type-2 AGN 

As for type-1 AGN, a relation between the intrinsic 15 ;um lu- 
minosity (L15) and the optical luminosity {Lr) for type-2 AGN 
is required to correct for the optical incompleteness. FLF04 has 



shown that such a relation exists and found it to be valid for the 
whole mid-IR population (excluding type- 1 sources) over more 
than 4 decades of mid-IR luminosity. In order to derive this re- 
lationship, optical identifications from RMS, ELAIS-S fields, 
HDF North (Aussel et al. 1999a, 1999b) and South (Mann et al. 
2002; Franceschini et al. 2003) were used. As opposed to type- 
1 sources, the ratio (Lis/Lr) has a dependence on L15 which 
can be expressed in the linear form 

log(Li5/LR) = 0.471ogLi5 - 5.02, 

where luminosities are expressed in solar units (see FLF04, 
Fig. 13). The data show an intrinsic Icr dispersion of 0.32. This 
equation implies that L15 increases ~3 times faster than Lr. The 
relation works for both nearby and distant objects (from z-0 up 
to z=0.7) and, therefore, has been used to estimate the fraction 
of sources lost due to the spectroscopic limit of the different 
surveys, i.e. the 0(z, L15) term used in the minimisation of the 
S function (see eq. 3). 

The LF was computed in the luminosity range logLis = 
[42,47] and in the redshift interval z = [0,0.7]. A total of 75 
sources (25 from ELAIS and 50 from RMS) were used to de- 
rive the LF shape and evolution. Best fit values for the LF pa- 
rameters, as well as the corresponding 2D-KS probabilities, can 
be found in Table 2 for each of the two representative mid-IR 
^-corrections assumed. 

For type-2 AGN, we followed the same fitting sequence 
as for type-1 AGN. The results for a fixed faint slope (PLE 
model) are found on Table 2 on rows labeled "C" (for the NGC 
1068 ^-correction) and "E" (for the Circinus ^-correction). In 
both cases, the value of the 2D-KS test (P2z)a:s=0.16 & 0.18 
respectively) does not reject the PLE model, but an important 
source overprediction is again present in the model if compared 
to the observed data for the low luminosity and moderate-z 
(z=[0. 1,0.6]) part of the sample. A better fit is found by assum- 
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Fig. 6. Observed and best fit LF for type-2 AGN at < z >-0.05 (continuous line) and < z >=0.35 (dash-dotted line). The solutions with a redshift dependent faint slope, asstiming a 
NGC 1068 ("D"; left panel) and Circinus ("F'; right panel) ^-corrections, are shown. Symbols and shaded areas as in Fig. 4. Due to the small redshift intervals used, and for the sake 
of clarity, the plots have been divided into two panels for the low-redshift (top) and intermediate-z sources (bottom). At intermediate-z, the thin line shows the low-z best fit. 



ing an evolving faint slope of the LF parameterised as specified 
in §3.2. In this case, in order to avoid unphysical solutions due 
to the low statistics at high redshift, the parameter was al- 
lowed to vary in the range [0, 6], while the parameter was al- 
lowed to vary in the range [0, 2]. Indeed, both parameters found 
a solution at the edges of the allowed ranges. The results can 
be found in Table 2 and Figure 6, fits labeled "D" and "F" for 
the NGC 1068 and Circinus ^-corrections respectively. 

Fig. 6 shows the observed and predicted space density dis- 
tribution of the type-2 AGN in two representative redshift in- 
tervals, low-z (z = [0,0.1]) where the IRAS sources dominate 
and intermediate-z (z = [0.1,0.6]), mainly populated by ISO 
sources. The overall best solution found for the two represen- 
tative SED, NGC 1068 (left panel, model "D") and Circinus 
(right panel, model "F"), are plotted. The observed space den- 
sity, errors, upper limits and the Icr confidence regions for the 
LF have been computed and plotted as described in §3.2. As 
observed for type-1 AGN, the errors introduced by the low 
statistics (light shaded area) are larger than the ones due to 
the assumptions on the Lis/Lr relation (dark shaded area). For 
example, in this case, while the error introduced by the as- 
sumed L\sILr relation in the determination of the density of 
type-2 AGN is -20% at the break luminosity (L\^) of the LF, 
the Umited statistics introduce, instead, a ~60% uncertainty. 
Therefore, also for type-2 AGN, in order to compute the er- 
rors on the LF parameters and on any derived quantity (counts 
and redshift distributions) the uncertainties introduced by the 
assumptions on the value of the Lh/Lr relation have been ne- 
glected. 

The rate of evolution found for type-2 AGN depends 
strongly on the assumed SED. If an NGC 1068 /r-correction is 
adopted, the rate of evolution is ~ 2.0 ± 0.5, lower than 



the one measured for type-1 sources (ki ~ 2.8 ± 0.3). On 
the other hand, an evolution rate ki ~ 2.5 + 0.5 is obtained 
if the Circinus SED is assumed, closer to the one measured 
for type-1 AGN. These differences in ki are caused by the 
stronger ^-correction introduced by the Circinus SED (Fig. 1), 
which generates higher intrinsic luminosities in the redshift 
range [0.1,0.6] compared to those given by the NGC 1068 k- 
correction. At variance, the impact of the assumed SED on all 
the other fitted parameters is small (see Table 2 & Fig. 6). 

For each best fit LF of the type-2 AGN, the predicted counts 
and redshift distributions were derived and compared with the 
observed data. Predictions at 15 jum for the integral counts are 
given in the left panel of Fig. 7, while the derived redshift 
distribution is shown in the right panel of the same figure. 
In both panels of Fig. 7 a thick-dashed line represents the 
Circinus fixed slope model ("E"), a thick-continuous line the 
Circinus variable slope model ("F"), a thin-dashed line the 
NGC 1068 fixed slope model ("C") and a thick-continuous fine 
the NGC 1068 variable slope model ("D"). 

A good agreement, within the errors, between the models 
and the data is observed. At variance with the solutions found 
for the type-1 AGN (§3.2), although a higher 2D-KS probabil- 
ity is obtained using a variable faint slope model, no major dif- 
ferences are observed in the integral counts down the faintest 
ELAIS flux (~lmJy; log5i5 ~ 0). The differences between 
the various models only become significant at fainter fluxes 
(5 15 < 1 mJy), especially at the fluxes reached by the Spitzer 
space observatory (5 24 ~ 0.01 mJy; see §6). Indeed, the rel- 
ative errors are rather uniform, with mean values around 30- 
50%. Therefore, the higher observed density of the deep ISO 
data (a factor ~10) is statistical significant. However, as previ- 
ously mentioned (§2), the higher density of AGN found in these 
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Fig. 7. Left ) Observed and predicted integral counts for type-2 sources using the the two representative SED: NGC 1 068 ii\im-dashed for the ''C" model and iMm-conlinuous for the 
''D" model) and Circinus {ihic]^-dashed for the "E" model and \h\c]^-conlimtous for the "F" model). Right) Observed and predicted type-2 AGN redshift distribution for the different 
models mentioned before. Lines as in left panel. Symbols and shaded area in both panels as in Fig. 5. Errors reported in the lower section of both panels have been computed and 
plotted as described in Fig, 5. 



fields is partly due to the classification method used, which is 
not only based on a pure optical spectral classification. The pre- 
dicted redshift distributions agree, within the errors, in the red- 
shift interval where the LF was computed {z - [0,0.7]). The 
difl'erent behaviour observed at z > 0.7 is understood by the 
different ^-corrections obtained from the adopted SED. The es- 
timated errors range from a factor ~ 1.3 at low redshift to ~2.0 
at the Umiting redshift used to fit the LF, and increase drasti- 
cally when the predictions are extrapolated at higher redshifts. 

3.4. The missing fraction of AGN 

In the previous sections we have seen that a flattening with red- 
shift of the faint slope of both type-1 and type-2 luminosity 
functions is slightly favoured by the data. This flattening would 
disappear if ~7 additional type-1 AGN, in the luminosity- 
redshift range logLi5=[42,45] and z=[0.2,3.2], and ~6 addi- 
tional type-2 AGN within logLi5= [42,44] and z=[0. 1,0.7] were 
found. We briefly discuss now two possible origins of incom- 
pleteness for the mid-lR AGN sample and their eff'ects on 
the observed LF: i) objects not identified due to the spectro- 
photometric limit of completeness of the survey and, ii) possi- 
ble tnisclassification of the objects. 

The spectroscopic incompleteness at faint optical mag- 
nitudes of the ELAIS sample is represented by the factor 
0(z,Li5) (§3.1). The corrections introduced through this term 
are represented in Fig. 8 as a function of redshift and mid- 
IR luminosity. The shaded areas in the figure outline the re- 
gion where type-1 (top) and type-2 sources (bottom), fainter 
than the spectro-photometric limit of the surveys, are expected 
to be found. A darker area in the plot indicates that a higher 
fraction of correction has been applied (i.e. higher number of 
lost sources) and it spans from ~1% (light-grey area) to ~30% 
(darkest area) of the total sky area available in a given inter- 
val of redshift and luminosity. The correction computed corre- 
sponds to ~4 type- 1 and ~ 1 type-2 sources in the L15 -z interval 



used in the minimisation process. If the suggested flattening is 
really due to an incompleteness effect, then the possibly missed 
objects are expected to follow an optical-mid-lR relation differ- 
ent from the one measured for the identified fraction of the high 
redshift 750 sources and in the local universe by //MS (z < 0.1) 
and adopted in our estimate of 0(z, iis). In the case of type-1 
AGN, if the average value of the assumed log(Li5/Ls) relation 
is increased by 2.5 times the observed dispersion of the relation 
(i.e. ~0.86 instead of 0.23), ~7 additional sources would be ex- 
pected at low luminosity and high redshifts, thus avoiding the 
need for a flattening of the LF. This higher Lis/Lr value could 
correspond, for example, to sources with a significantly larger 
amount of gas and dust with respect to the local samples of 
Elvis et al. (1994) and Spinoglio et al. (1995). For type-2 AGN, 
the relation derived by La Franca et al. (2004) implies that for 
the low luminosity (L15 < L^j) and moderate redshift (z ~0. 1- 
0.5) sources, all the possible optical counterparts should have 
been observed. Only the high luminosity (L15 > ^15), higher 
redshift (z > 0.5) sources are aff'ected by incompleteness. Also 
in this case a significant higher value of the Lis/Lr relation 
would be required in order to avoid the need for a flattening of 
the LF 

A second source of incompleteness could be caused by a 
spectroscopic mis-classification of the sources. The increase in 
the mean redshift of the ISO population with respect to the lo- 
cal IRAS sample (see Table 1) implies that a larger fraction of 
the host galaxy light is collected for a given size of the sUt 
used for the optical spectroscopic identifications. The conse- 
quence is a dilution of AGN signatures that aff'ects the mea- 
sured equivalent widths (EW) and line ratios (see e.g. Moran 
et al. 2002). The effect is less important for type-1 sources due 
to the presence of easily identified broad emission lines caused 
by a direct view of the central energy source with a low level 
of obscuration. The situation is more complicated for type-2 
AGN since their higher degree of obscuration suppresses an 
important fraction of the optical light from the nuclei. The host 
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Fig. 8. Luminosity-redshift distribution of tiie type-i (Uip) and type-2 (hoiiciu) AGN. 
Shaded areas represents the fraction of incompleteness due to the optical spectroscopic 
limits. Contours range from 75% (black) to 95% (light-gray) of the total area lost and 
have a 5% step between them. For type-1 AGN ~4 missed sources are expected within 
Z=[0.1,1.2]. IVIost missed type-2 AGN are only expected with logLis > 44.5 nt z > 0.7, 
a redshift range outside the current analysis (see text). Only ~ 1 missed type-2 source is 
expected to lie within z = [0.3,0.5] and Z.15 = [43.5,45.0]. Dotted lines in both panels 
show three representative flux limits at 1, 4 and 300 mjy. 



galaxy would then naturally contribute with a larger fraction to 

the observed optical light from the nuclear source, sometimes 
even dominating it. The effect is maximised for the intrinsi- 
cally faint AGN, whose optical features will be highly diluted 
into the host galaxy spectrum. This effect could be relevant in 
producing the observed flattening at high redshift of the lumi- 
nosity function. 

Observations with XMM and Chandra have uncovered a not 
negligible number of AGN for which the classification given by 
their optical spectra does not correspond to the X-ray classifi- 
cation of the source (e.g. Fiore et al. 2000, 2003; Barger et al. 
2005). In their analysis of deep (~70 Ksec long) XMM observa- 
tions on ELAIS SI -5, La Franca et al. (in preparation) find that 
the fraction of misclassified type-1 AGN appears to be negli- 
gible. Vice-versa, in the case of type-2 AGN, the observations 
with XMM have revealed that about 1 0-20% of the ISO counter- 
parts optically classified as no-AGN (i.e. normal and starburst 
galaxies) show X-ray luminosity consistent with AGN activ- 
ity. This fraction, in the case of ELAIS-S, would imply that as 



Luminosities have been converted to 15 fim using the adopted SEDs. 
W From Xuetal. 2003. 

Best tit solutions for type-1 (model "B") and type-2 (Circintis SED, model "F"). 
Intermediate redshift, z=0.6. 



many as 10 to 25 sources, classified in the optical as starburst 
or normal galaxies, may harbour an AGN that can contribute to 
the observed mid-IR flux. If this is the case, the number of true 
type-2 AGN may be higher by a factor ~ 1 .4 - 2.0 than that used 
in this analysis and this would have important consequences for 
the observed LF and its evolution. 

We conclude that the observed flattening with redshift of 
the LF for low luminosity AGN can be explained if we assume 
that the incompleteness in the ELAIS sample was not properly 
modeled, for two possible reasons: (i) the unidentified fraction 
of sources may have a L\s/Lr ratio different from that of the 
identified AGN or, (ii) more likely, a relevant fraction of type- 
2 AGN could have been mis-classified due to dilution of the 
optical nuclear spectra by the hosting galaxy. These sources of 
incompleteness are probably almost negligible for type-1 AGN 
(especially the last one), while our measure of the density of 
type-2 AGN has instead to be considered a lower limit. 



4. Comparison with previous works 

4.1. The IRAS 12nm local luminosity function 

A local luminosity function (LLF) for type-1 and type-2 AGN 
was derived at 12/im by RMS. Their results can be directly 
compared to our findings at z=0. The values for the LF pa- 
rameters reported by RMS and those computed by us are sum- 
marised in Table 3. 

Two parameters of the LF computed by RMS, the faint 
slope (a) and the luminosity break (L'^^), difi"er significantly 
from our solutions for both the type-1 and the type-2 AGN. 
The difference on the value of the luminosity break can not be 
explained by the SEDs conversion factors from 12 to 15 fim, 
since they imply a small change in the intrinsic luminosity of 
the sources (A logL[12 - 15/im] ~0.0 for type-1 and ~0.1 for 
type-2 AGN). These discrepancies can instead be understood 
if we take into account that: (i) the faint slope was fixed a pri- 
ori in the fitting procedure by RMS; (ii) the strong effect of 
evolution in the different luminosity bins was not considered 
by RMS and, (iii) our best fit is obtained by taking into ac- 
count not only the local IRAS sources but also the high redshift 
ISOCAM population. 
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4.2. The IRAS 25nm combined AGN sample 

A local combined (i.e. type-1 + type-2) AGN LF has been de- 
rived by Xu et al. (2003) at 25 fim. It is based on an IRAS colour 
selected (S eo/jm / S 25i.im < 5) sample obtained from a complete 
catalog compiled by Shupe et al. (1998). Assuming a luminos- 
ity evolution equal to the optical one (Boyle et al. 2000) and a 
set of SEDs with a dependence on luminosity, Xu et al. (2003) 
were able to derive the predictions for the AGN counts, the red- 
shift distribution and the contribution to the cosmic background 
from the UV to the sub-mm wavelengths. A direct comparison 
of the two LFs can then be made by translating our 15 fim lu- 
minosities into 25 /xm luminosities, using the adopted SEDs for 
each source class. The results of this comparison, in two differ- 
ent redshift intervals, can be found in Fig. 9 and Table 3. We 
find a good agreement between the two LFs (dashed line for 
the 25//m selected sample; continuous line for our type-1 + 
type-2 15 jum LF) in the local universe (z~0.1). The difference 
found at intermediate-z (z=0.6, the maximum z used to derive 
the type-2 LF) and low luminosities L < L*) is mainly due to 
the observed decUne in space density at this redshift for the ob- 
scured sources. The difference at higher luminosities (L > L*) 
can be explained by the stronger evolutionary rate found by 
Boyle et al. (2000) for optical selected type-1 AGN (ki ~ 3.5) 
and used by Xu et al. (2003) in their analysis with respect to 
the values found in this work, k ~ 2.0 - 2.9 depending on the 
AGN type. The difference in the LF induced by these different 
kt parameters is already significant at z ~ 0.6. 

4.3. ISOCAM Deep observations 

The deep observations carried out by the ISOCAM instrument 
on-board of ISO place interesting constraints on the shape and 
behaviour of the LF. For type- 1 AGN they favour a PLE model 
scenario (dashed line in Fig. 5, left panel). This suggests that 
a small fraction of the low-luminosity ELAIS-S type-1 AGN 
may still be hidden within the rest of the mid-lR population. If 
this is the case, they do not follow the optical-mid-lR luminos- 
ity relation assumed for the rest of the type-1 sources as already 
discussed in §3.4. 

The observations of type-2 sources in these fields indicate 
that all our best fit solutions underestimate them by more than 
1.5£r. This result is not surprising since most of the type-2 
sources identified in these surveys have been classified using 
criteria different from the pure optical spectral classification 
(i.e. SED reconstruction, X-ray properties and radio emission). 
Indeed, this was the reason for which we did not consider these 
sources in our minimisation process. A comparison between 
their observed space density and our predictions has to be made 
with care due to the different classification methods used. 

We would also like to note here that even if cosmic variance 
can increase the HDF-N, -S and CFRS quoted errors in Fig. 
5 and Fig. 7 (by as much as 20-100% in such a small fields, 
Somerville et al. 2004), it can not fully justify the observed 
difference since in all these fields the observed counts are found 
to be systematically higher than our predictions. 




-m l 

42 43 44 45 46 47 
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Fig. 9. Comparison of LFs at 25 //m. The derived 15jum LF for type-1 (model "B") and 
type-2 (Circinus SED, model "F') AGN, converted to 25 yum luminosities, are compared 
with the global AGN LF fix>m Xu et al.(2003). The comparison is made in two different 
redshift intervals: z=0.1 (top panel) and z=0.6 (bottom panel). The values for the parame- 
ters of the LFs are reported in Table 3. 

5. Contribution to tlie CIRB 

The contribution of AGN to the intensity of the Cosmic 
Background Ught in the Infrared (CIRB) has been derived for 
all the models presented in Table 2. The intensity of the CIRB 
at 15 /zm for a given population has been computed as, 

I=^J dLisJ dz^5,(Li5,z)0(Li5,z) 

where Sv(Lis,z) represents the observed flux of a source with 
an intrinsic luminosity L15 at redshift z. The integration has 
been performed for log Li5(z=0)= [42,47] up to z=4. For type- 
2 sources the same redshift cutoff as for type-1 (ze«r=2.0) has 
been assumed. Column X in Table 2 summarises the expected 
background light for each model and the associated errors 
given by the Icr dispersion on the LF parameters as described 
in §3.2. 

The integrated light from type-1 AGN galaxies provides 
a contribution of (4.2-12. l)x 10"" Wm^^sr"' (in units of 
v/y) at 15jum. This value corresponds to (1.6^.5)% of the 
15/im background as reported by Elbaz et al. (2002) and 
confirmed by Metcalfe et al. (2003; 2.7 xlO"'' Wm^^sr"') 
using deep lensed observations with ISOCAM. Our predic- 
tions are in agreement with what previously estimated by M02 
(5.7 X 10-" Wm--sr-i). 

The contribution of type-2 sources is found to be (5.5- 
14.6)xl0"'^ Wm"^sr"\ depending significantly on the as- 
sumed SED and the parameterisation for the evolution of the 
faint slope [a(z)], and accounts for (2.0-5.4)% of the CIRB at 
15 fim. We note here that, even if the space density of the ob- 
scured sources, up to redshift 0.7 (where type-2 AGN are de- 
tected), is a factor 2-6 higher than that of the un-obscured ones 
(in the luminosity range L15 ~[43.5,45.0], around L*, where 
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most of the CIRB is produced), their contributions to the CIRB 
at 15 fim are very similar. The main reason is due to the fact that 
type-2 AGN have a strong ^-correction imposed by the mid-IR 
SED. Moreover, the effect is increased by a sUghtly lower evo- 
lution rate in comparison to the type-1 AGN, and a significant 
decline of the faint space density at higher z. 

The total estimated contribution of AGN represents ~4- 
10% of the total light observed in the mid-IR. This fraction 
is about half of the one ascribed to the AGN by M02. M02 
study was limited to type-1 AGN, and assumptions based on 
the unified model and the local ratio of type-2 to type-1 (~4; 
Maiolino & Rieke 1995) were made to derive their contribu- 
tion. The smaller ratio of type-2 to type-1 found in our fields 
can explain the observed difference in our estimates. 

Mid-IR studies based on IRAS and ISO observations have 
indirectly estimated a contribution of AGN to the CIRB not 
larger than 5-10% (e.g. Franceschini et al. 2001; Xu et al 2001, 
2003). This is the maximum room 'left' in their models by 
the strongly evolving starburst population. The total contri- 
bution was in any case uncertain since mid-IR selected type- 
2 AGN and starburst galaxies were treated as a single popu- 
lation. Another estimate of this contribution comes from the 
X-ray band (0.5-10 keV), which offers a better wavelength 
regime where to select and identify obscured sources (unless 
the sources are Compton-thick, Nh > 10^^ cm"-^). A cross- 
correlation of X-ray and IR sources detected by deep obser- 
vations in the Chandra Deep field North (CDFN; Brandt et 
al. 2001) and in the Lockman Hole (Hasinger et al. 2001) al- 
lowed Fadda et al. (2002) to estimate the maximum fraction of 
the CIRB produced by AGN. This fraction is found to be of 
(17±6)%, a factor 2-4 higher than our predictions. 

On the basis of the hard X-ray determination of the LF 
and semi-empirical SEDs (linking the X-ray to the Infrared), 
Silva et al. (2004) have derived the contribution of AGN and 
their host galaxy to the CIRB. For type-1 AGN, the agreement 
between their results (~ 7 x 10"^' Wm^^sr ') and the ones 
presented here is very good. 

There is instead a significant difference for type-2 sources. 
The predictions from Silva et al. (2004; 3 x 10"^° Wm ^sr ^) 
are, in the most favourable case (model "E"), at least a factor 
of 2 higher than ours and closer to the numbers provided by 
Fadda et al. (2002). The higher efficiency of the hard X-ray 
observations to select obscured AGN can explain the discrep- 
ancy between our results and the ones presented by Fadda et al. 
(2002) and Silva et al. (2004). 

The conclusion drawn from the comparison with the above 
mentioned studies is that, as discussed in §3.4, it is likely that a 
significant fraction of type-2 AGN is hidden within the rest of 
the mid-IR selected population, the normal and starburst galax- 
ies. 

6. Mid-IR counts predictions for Spitzer 

The Spitzer space telescope is currently performing several 
deep mid- and far-infrared observations in selected areas of the 
sky. Following the results presented in this paper we have esti- 
mated the expected number of AGN sources as a function of the 
flux in the Spitzer mid-IR band. The derived integral counts for 



the MIPS instrument at 24 /um are shown in Figure 10. The re- 
sults are given in three panels that correspond to three different 
combinations of type-1 (thick dashed) and type-2 (thick dot- 
dashed) models (Table 2). The total AGN contribution (type- 
1 + type-2) is given by a thick continuous line. Our estimates 
are compared to the total contribution of AGN from the mod- 
els of Xu et al. (2003) and Silva et al. (2004). These models 
are plotted as a thin triple-dot-dashed (Xu et al. 2003) and thin 
continuous (Silva et al. 2004) lines. In the case of the Silva et 
al. (2004) model the individual contributions from type-1 (thin 
dashed line) and type-2 (thin dotted-dashedline) AGN are also 
shown. Best fit solutions with a redshift dependence on the faint 
slope (fits "B" and "F" for type-1 and type-2 respectively) are 
plotted in the left panel. The central panel presents the expecta- 
tions assuming the best fits without any evolution for the faint 
slope of the LF (PLE models, fits "A" and "E"). Finally, the 
combination of type-1 and type-2 models that better fits the ex- 
pectations from Silva et al. 2004 (fits "B" for type-1 and "E" 
for type-2) is given in the right panel. 

At a 24yum flux of 0.01 mJy the total AGN counts given 
by our best fit solutions underestimate the predictions by Silva 
et al. (2004) and Xu et al. (2003) by a factor of -2-3 (Fig. 10, 
left). With these models we predict to find ~1200^3qq type-1 
and ~1000;|;25Q type-2 optically classified AGN per sq. degree 
down to a Spitzer flux Umit of 524^m = 0.01 mJy. If compared 
to the individual contributions of type-1 and type-2 AGN given 
by the Silva et al. (2004) model, we note that the large disagree- 
ment is mainly due to the low number of type-2 AGN expected 
in our models (a factor ~5 lower) and caused by the rapid flat- 
tening with z of the faint LF slope. The large effect of the LF 
faint slope decUne in the counts is evident if compared to the 
PLE models predictions plotted in the central panel. The PLE 
models overestimate the 24 jim counts by both Xu et al. (2003) 
and Silva et al. (2004) below 0.06 mJy. A combination of the 
model "B" (variable LF faint slope for type-1 AGN) and model 
"E" (fixed faint slope model for type-2 AGN) produces the best 
approximation to the counts derived by the Silva et al. (2004) 
model based on the AGN hard X-ray LF (Fig. 10, right). If this 
is the case, then ~1200 type-1 and ~5600 type-2 AGN per sq. 
degree are expected brighter than 0.01 mJy at 24 pm. 

We want to note here that: i) although the expected 
Spitzer/24pm integral counts derived by our different models 
differ by factors as large as ~6-10, they are very similar down 
to the flux limits of the ELAIS-South Surveys (shaded area in 
Fig. 10), responsible for the majority of the sources observed at 
high redshift in our sample, ii) even if some of our models pro- 
duce counts comparable to the ones derived from X-ray LFs at 
SiA/um = 0.01 mJy (Fig. 10, central and right panel), they are 
always lower at brighter fluxes. This second point explains the 
observed difference (by a factor of ~3) between our models and 
Silva et al. (2004) in the contribution of AGN to the CIRB. 

7. Conclusions 

Combining ISO/15 pm observations in ELAIS flelds, HDF-N, 
HDF-S and the CFRS with the local IR population detected by 
IRAS at 12 pm, we have derived the evolution for the mid-IR 
selected and optically classified AGN. While a similar study 
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Fig- 10. The expected Spilzer MIPS integral counts for AGN at 24 /im based on our best fit models (thick-dashed for type-1, thick-detl-dashed for type-2 and thick-conlinuems line 
for tile total). Tiie panels represent dilTerent combinations of type-1 and type-2 models from Table 2. Left Models ''B" and "F" for type-1 and type-2 respectively. Central Models "A" 
and "E". Right Models "B" and "E". For comparison, the expectations from the dilferent AGN classes derived by Silva et al. (2004) are also shown (Type-1: thin dashed line; Type-2: 
thin dot-dashed; Total: thin continuous). The total AGN contribution by Xu et al. (2003) is represented by a thin triple-dot-dashed line. The fight-grey shaded areas represent the 
expected flux range at 24j:(m of a source with a ISpim flux of 1 mjy given the assumed type-1 and type-2 SED. 



has already been done for QSOs-i-Seyfert-1 sources in the past 
(Matute et al. 2002), we have presented here for the first time 
the rate of evolution shown by the obscured, type-2 AGN. Our 
results are briefly summarised as: 

- Type-1 AGN evolve following a double power law LF and 
a rate of evolution of the form (l+z)'"- with kL ~2.9. The 
data are consistent, within the errors, with a PLE model. 
However, a flattening at high redshift of the faint luminosity 
end of the luminosity function is marginally favoured by the 
data. 

- Type-2 sources evolve with a slightly lower rate than type- 
1, with^i ranging from ~1.8 to ~2.6 depending strongly on 
the assumed SED used to compute the K-correction. The 
best fit solution favours a luminosity dependent luminosity 
evolution (LDLE) model, while the PLE model is statisti- 
cally acceptable. However, integral counts reported for the 
deep 15 /zm observations in the HDF-N, -S and CFRS sug- 
gest an evolutionary scenario closer to the PLE model. We 
expect a significant number of type-2 AGN to be hidden 
within the optically classified normal and starburst galax- 
ies. 

- In the volume of the universe where both types of sources 
(type-1 & type-2) are observed (z=[0,0.7]) the ratio of type- 
2 to type-1 is ~2-6. This value is in agreement, within the 
errors, with the one measured from the local optically se- 
lected samples, where a ratio close to 4 is found (Maiolino 
& Rieke 1995). 

- The total contribution of AGN to the CIRB at 15;um is of 
the order of ~4-10% of the total background light mea- 
sured by Metcalfe et al. (2003), divided approximately 
equally between the two classes (type-1 & type-2). The 
contribution of type-1 sources is in good agreement with 
the results presented in previous studies (Matute et al. 



2002). Comparison with results from X-ray selected sam- 
ples (e.g. Fadda et al. 2002, Silva et al. 2004) shows fliat a 
significant fraction of obscured type-2 sources are missed 
in our sample. We argued that the optical spectral classifi- 
cation of the mid-IR sample, and not the mid-IR selection, 
is the principal responsible for the missed fraction. 
- Estimates for the number of mid-IR selected and optically 
classified AGN, expected from Spitzer at 24 foa, are given 
according to the best model results. We expect ~1200 type- 
1 and ~1000 type-2 optically classified AGN per sq. de- 
gree at a flux limit of S24jjm = 0.01 mJy. Our 24 /zm counts 
derived for type-1 sources agree very well with previous 
works. At a flux limit of S24iim = 0.01 mJy the expected 
counts for the obscured population from our best fit mod- 
els are a factor ~5 lower than the expectations from models 
based on X-ray LFs (e.g Silva et al. 2004). 

The results presented here for type-1 sources are quite ro- 
bust since a high completeness is expected for these sources. 
Unlike type-1, type-2 AGN show a larger spread in their mid- 
IR and optical properties and a significant fraction can be mis- 
classified as starburst or normal galaxies. As a consequence, 
the results for type-2 AGN derived here can only be consid- 
ered as a lower limit to their true density and a first approxima- 
tion to the evolution of these sources in the mid-IR. The true 
mid-IR space density of obscured sources has to be determined 
combining their mid-IR properties and the optical classification 
with the information available at other wavelengths, especially 
in the X-rays. 
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